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INTRODUCTION  (STATEMENT  OF  WORK) 

MOTIVATION 

Renewed  scientific  interest  has  arisen  throughout  the  world  as  to  the  potential 
application  of  MHD  processes  for  advancement  of  flight.  Specific  to  this  are  the 
collaborative  research  works  going  on  between  the  US  and  IVTAN  in  the  field  of 
electromagnetic  (EM)  Advanced  Flow  and  Flight  Control  for  hypersonic  aircraft  and 
single-stage  earth-to-orbit  vehicles.  Among  the  areas  of  interest  is  that  of  utilization  of 
MHD  as  a  means  for  enhancing  the  speed  and  range  of  scramjets  through  a  concept 
known  as  MHD  energy  by-pass  [1-4].  Currently,  scramjet  operation  is  limited  to  free 
stream  flight  Mach  numbers  around  10.  This  limitation  arises  from  the  excessive 
temperature  that  develops  at  high  Mach  conditions  as  a  result  of  the  slowing  down  of 
the  propulsion  drafted  hypersonic  air  stream  within  the  propulsion  system  inlet.  High 
temperature  air  (>  1,600K)  entering  the  scramjet  combustor  inhibits  combustion  due  to 
the  natural  dissociation  of  chemical  species  that  occurs. 

The  MHD  energy-bypass  concept  is  a  two  stage  process.  First  it  employs  an  MHD 
generator  situated  in  the  inlet  to  extract  enthalpy  from  the  propulsion  air  stream  to  cool 
the  flow  prior  to  entering  the  scramjet.  This  extracted  thermal  energy  is  converted  to 
electrical  power  through  the  MHD  generator  process  of  electromagnetic  induction.  With 
the  scramjet  combustor  incoming  flow  cooled,  the  flight  corridor  of  the  aircraft  can  be 
enhanced  to  higher  Mach  number.  In  the  second  stage  of  this  process,  the  electric  power 
extracted  by  the  generator  is  used  to  power  an  MHD  accelerator  in  the  propulsion 
systems  tailpipe.  The  tailpipe  MHD  accelerator  is  the  generator  process  in  reverse 
(motor  mode  of  operation)  and  results  in  adding  energy  directly  back  to  the  stream  flow 
as  kinetic  energy  through  the  action  of  Lorentz  forces.  Consequently,  the  flow  energy 
extracted  in  the  inlet  generator  is  recovered  in  the  tailpipe  MHD  accelerator, 
negating/cancelling  dynamic  stream  losses. 

The  figure  1  exemplifies  the  MHD  energy-by  pass  concept. 


Flow  Energy 


4  5 


MHD  Scramjet  MHD 
Generator  Accelerator 


Fig-1.  Schematic  of  MHD  energy  by  pass  concept. 
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A  potential  area  of  collaboration  has  recently  arisen  between  IVTAN  and  NASA’s 
Advanced  Space  Transportation  Program  (ASTP).  NASA  ASTP  is  interested  in  this 
propulsion  concept  as  a  potential  new  system  for  single-stage-to-orbit  vehicles. 
Collaborative  studies  done  by  Dr.  V.A.Bityurin  of  IVTAN  and  a  U.S.  Small  Company, 
LyTec,  LLC  have  screened  the  thermodynamics  of  the  MHD  energy  by-pass  concept 
[5-7].  The  promising  results  of  that  work  have  led  to  increased  interest  by  NASA  ASTP 
to  explore  operation  of  components  of  this  system  in  experiments. 

One  of  the  key,  and  least  understood  components  of  the  system  is  the  MHD  accelerator. 
Expertise  and  experience  in  this  area  does  exist  at  IVTAN  TsAGI  and  our  US 
counterparts  recognize  our  unique  knowledge  in  this  field.  In  addition,  IVTAN  has 
readily  available  access  to  an  operational  MHD  accelerator  that  can  be  used  as  a  test 
bed.  In  anticipation  of  the  follow-up  collaborative  program  with  NASA  ASTP 
(currently  being  requested  through  the  agency’s  international  programs),  we  have 
proposed  to  EOARD  an  initial  effort  in  the  exploration  of  MHD  acceleration  for 
aerospace.  Our  proposed  accelerator  project  is  equally  applicable  to  exploring  MHD 
propulsion  and  to  utilizing  MHD  acceleration  as  a  flow  driver  for  hypersonic  ground 
testing. 

Our  intention  under  this  proposed  project  was  to  conduct  some  initial  feasibility  and 
analysis  studies  on  application  of  the  MHD  accelerator  for  propulsion  and  flow 
acceleration.  In  this  analytical  effort,  we  plan  to  analyze  the  MHD  acceleration  process 
as  a  means  to  improve  our  understanding  of  this  type  of  device  and  to  qualify  its 
application  to  MHD  energy  by-pass  propulsion.  This  work  is  a  prelude  to  the 
collaborative  IVTAN-NASA  follow-on  program  that  will  be  both  experimental  and 
theoretical.  In  addition,  we  use  this  proposed  EOARD  work  as  a  vehicle  to  make  the 
experimental  MHD  accelerator  facility  operational.  The  facility  we  haves  utilized  is  an 
MHD  accelerator  driven  wind  tunnel  that  been  idle  since  1995.  We  envisioned  that  it 
would  require  approximately  two  and  one  half  months  of  technician  and  engineering 
effort  to  refurbish  the  tunnel  to  operational  status.  This  anticipation  has  been  fully 
confirmed  by  the  proposed  work  done. 

Our  proposed  project  will  be  to  the  benefit  of  all  parties  involved,  present  and  future. 
The  MHD  accelerator  test  facility  note  operational  provides  a  working  prototype  that 
can  be  directly  researched  to  qualify  the  MHD  accelerator  for  potential  application  to 
advanced  propulsion  concepts  of  interest  to  NASA  and  others.  This  laboratory  readiness 
effort  makes  available  a  hypervelocity  wind  tunnel  that  can  be  used  as  a  test  bed  by 
others  (US  Air  Force  and  US  industries)  in  the  wider  range  of  programs  and  initiatives 
that  are  taking  place  in  exploration  of  EM  Advanced  Flow  and  Flight  Control.  Indeed, 
our  proposed  project  falls  well  within  the  mission  of  the  EOARD  of  expanding  the  U.S. 
and  Russian  scientific  cooperation  and  resources  for  revolutionary,  breakthrough 
aeronautical  research. 
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PROJECT  OBJECTIVES 

The  overall  objective  of  our  proposed  project  was  to  make  ready  (“wet”)  an 
experimental  MHD  accelerator  laboratory  that  can  be  used  in  the  future  for  both 
experimentation  into  MHD  propulsion  and  as  a  hypervelocity  wind  tunnel  for 
experiments  into  hypersonic  flow  and  flight  control  systems.  Within  the  confines  of  our 
proposed  program  we  seek  to  achieve  the  following  specific  technical  objectives: 

1.  To  analysis  and  report  on  the  potential/capability  of  MHD  acceleration  assisted 
advanced  propulsion  and  to  provide  quasi-steady  state  hypervelocity  ground  tests; 

2.  To  ready  an  existing  experimental  MHD  accelerator  laboratory  for  experimental 
evaluations; 

3.  To  develop  an  experimental  plan  for  conduct  of  an  open  demonstration  experiment. 
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PROJECT  PLAN 

The  following  Table  provides  the  work  breakdown  structure  and  brief  descriptions  of 
the  discrete  tasks  of  our  Project  Plan. 


1.0  Background  Studies 

•  Conduct  Literature  Reviews  on  the  MHD 

Acceleration  Process  and  Past  Experiments 

2.0  Analysis 

•  Conduct  Literature  Reviews  on  MHD  acceleration 
and  Past  Experimentation 

•  Develop  Methods  for  Analysis  of  MHD  Accelerators 

•  Provide  Analytical  Support  to  Experimental 

Planning  and  Interpretation  of  Experimental  Results 

3.0  Lab  Qualification  and 
Preparation 

•  Qualify  Test  Worthiness  of  Available  Laboratory 
Hardware  and  Equipment 

•  Re-furbish  Laboratory  to  Operational  Status 

•  Define  and  Implement  Diagnostics  and  Support  Data 
Systems 

•  Define,  Design,  Fabricate  Test  Hardware  Models 

4.0  Experimental  Plan 

•  Define  Basic  Experimental  Approach(s)  -  Test 
Conditions  and  Run  Scenarios 

•  Develop  Pre-Test  Analytical  Predictions 

•  Develop  a  Detailed  Test  Plan  for  Conduct  of  a  Open 
Demonstration  Experiment 

5.0  Program  Plan 

•  Develop  an  Overall  Program  Plan  for  Conduct  of 
Demonstration  Experiment 

•  Co-ordinate  Program  Plan  with  Potential  Facility 

Users 

6.0  Reporting 

•  Draft  a  Final  Report  Detailing  the  Facility  and  Its 
Capabilities  as  Project  Deliverable 
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FACILITIES  AND  EQUIPMENT 

A  unique  MHD  accelerator  of  TsAGI  laboratory  headed  by  prof.  V.  Alfyorov  has  been 
made  available  through  IVTAN.  The  laboratory  has  been  made  operational  within  a  few 
months  after  project  initiation  through  standard  re-furbishment  work  under  Task  3.0. 
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I.  CURRENT  STATE  AND  POTENTIALITIES  OF  HYPERSONIC 
MHD-GAS  ACCELERATION  WIND  TUNNELS 


1.  INTRODUCTION 

Developments  of  hypersonic  wind  tunnels  with  MHD-gas  accelerators  were  initiated 
late  in  the  1950s  when  limited  capabilities  of  classical  hypersonic  wind  tunnels  based  on 
the  test  gas  acceleration  due  to  its  adiabatic  expansion  became  obvious.  These  studies 
were  started  almost  at  the  same  time  in  the  USA  (NASA,  Langley  Research  Center, 
AEDS,  General  Electric  Company)  and  in  Russia  (TsAGI).  Some  works  were  carried 
out  in  France  (ONERA)  and  in  Germany. 

The  present  review  analyzes  different  designs  of  such  facilities.  Main  characteristics  of 
the  hypersonic  wind  tunnel  of  TsAGI  and  the  results  of  its  application  to  solve  practical 
problems  of  aerodynamics  are  considered.  The  flow  conditions  over  bodies  in  a  free 
flight  are  compared  with  those  realized  experimentally  in  the  test  section  of  the  facility 
at  the  same  flow  velocities  of  F~8km/s.  The  advantage  of  applying  MHD-gas 
acceleration  facilities  to  solve  problems  of  developing  and  testing  scramjets  for  TAV 
vehicles  is  demonstrated  and  the  data  on  possible  gas  dynamic  and  electrodynamic 
parameters  of  these  facilities  and  their  design  versions  are  presented.  Finally,  basic 
problems  to  be  solved  in  developing  hypersonic  MHD-gas  acceleration  facilities  are 
listed. 

The  aerospace  engineering  was  progressing  simul-taneously  with  the  development  of 
facilities  intended  for  testing  novel  vehicles  or  their  models  in  ground  laboratory  stands. 

Before  the  middle  50s,  hypersonic  wind  tunnels  were  based  on  the  principle  of 
converting  the  total  stagnation  enthalpy  into  the  kinetic  flow  energy  during  gas 
expansion  in  a  supersonic  nozzle.  The  programs  of  piloted  Earth  reentry  after  the  flight 
to  the  Moon  in  early  60s  have  posed  the  problem  of  constructing  facilities  capable  of 
reproducing  flow  parameters  at  M-20-25  and  Re~107-108  with  an  energy  flow  density 
of  ~107m2/s2.  At  that  time,  it  became  clear  that  the  capabilities  of  classical  wind  tunnels 
to  solve  this  problem  were  limited.  An  important  limitation  on  obtaining  high 
parameters  of  the  flow  with  a  velocity  of  F~5-7km/s  was  the  attainment  of  stagnation 
parameters  Po>108Pa  and  To>l 0000-1 2000K  in  the  wind  tunnel  plenum  chamber  when 
growth  of  radiation  losses  (gas  radiates  as  a  black  body)  and  break-down  of  the  nozzle 
throat  for  the  time  less  than  10~3s  are  governing  factors.  Besides,  high  stagnation 
parameters  led  inevitably  to  non-equilibrium  gas  flows  in  hypersonic  nozzles. 

A  way  of  increasing  flight  parameters  to  be  simulated  without  achieving  ultrahigh 
parameters  in  the  facility  plenum  chamber  is  the  method  of  magnetogasdynamic  gas 
acceleration  enabling  a  direct  increase  in  the  kinetic  flow  energy  under  the  action  of  the 
Lorentz  forces. Various  MHD-accelerators  were  also  proposed  for  wind  tunnels,  among 
which  were  accelerators  with  a  travelling  electromagnetic  wave,  Hall-type  and  Faraday- 
type  accelerators.  Each  of  them  was  characterized  by  its  own  operating  range  of  flow 
density.  For  the  first  accelerator,  it  was  p=10”7-10“4kg/m3,  for  the  second  one  p=10'6- 
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10'2kg/m3,  and  finally,  for  the  Faraday-type  one  p=10_2-10kg/m3.  The  Faraday-type 
MHD-accelerators  implying  the  Hall  current  neutralization  by  means  of  electrode 
sectionalization  received  the  widest  recognition  because  the  facilities  providing 
maximum  possible  high  flow  densities  are  of  particular  interest  for  aerodynamic 
investigations.  The  simulation  of  hypersonic  flights  necessitates  a  correct  reproduction 
of  the  following  conditions: 


1 .  Flight  velocity 

2.  Mach  number 

3.  Small  initial  oxygen  dissociation  level 

4.  Correct  Damkohler  number  correlation  for  flows  in  the  test  section  and  in  full-scale 
conditions 

5.  Rather  great  run  duration  to  solve  problems  of  experimental  investigations. 

The  choice  of  basic  accelerator  characteristics  depends  on  test  objectives. 


Essentially  all  hypersonic  MHD-gas  acceleration  wind  tunnels  being  developed  in  the 
early  60s  have  the  following  conceptional  design  (Fig.2). 


Fig.2.  Schematic  of  MHD-gas  acceleration  wind  tunnel.  1  -  arc-heater,  2  -  mixing 
chamber,  3  -  metering  device,  4  -  primary  supersonic  nozzle,  5  -  MHD-accelerator, 

6  -  secondary  nozzle,  7  -  test  section 


A  gas  (air,  nitrogen  etc.)  heated  in  the  source  of  a  hot  gas  up  to  7o~3500-4000K  and 
added  with  alkali  seeds  expands  in  a  primary  supersonic  nozzle  up  to  M=l. 8-2.5  to 
enter  a  sectionalized  Faraday-type  MHD-channel.  The  quantity  of  electrodes  with 
independent  power  sources  amounts  to  -30-120.  The  magnetic  field  with  induction  of 
5-1 .5-5T  is  applied  perpendicularly  to  the  flow  so  that  to  provide  the  Lorentz  force 
accelerating  the  gas.  In  a  secondary  hypersonic  nozzle  behind  the  MHD-channel,  the 
Mach  number  increases  from  3-4  to  15-20  as  a  result  of  adiabatic  cooling  and 
conversion  of  static  enthalpy  into  velocity.  Test  models  are  installed  in  the  test  section, 
while  supersonic  and  subsonic  diffiisors,  a  cooler  and  an  exhaust  system  are  mounted 
behind  the  test  section.  Either  a  vacuum  reservoir  or  an  ejector  system  are  used  as  an 
exhaust  system.  Fig.2  shows  the  schematic  of  the  facility. 

Below  consideration  is  given  to  physical  processes  proceeding  in  the  components  of  the 
above  outlined  facility  and  to  designs  of  its  units  used  in  different  organizations. 


13 


Final  Report  to  EOARD  from  IVTAN 
10.07.2000 


2.  MAIN  COMPONENTS  OF  AN  MHD  WIND  TUNNEL 

2.1.  Hot  Gas  Source 

A  hot  gas  source  should  provide  required  temperature  and  pressure  levels  sufficient  to 
achieve  a  high-temperature  flow  with  a  minimum  initial  gas  dissociation  degree  at  the 
primary  nozzle  exit.  Besides,  flow  fluctuations  in  it  must  be  minimized.  The  level  of 
stagnation  parameters  is  governed  by  a  characteristic  operation  time  and  a  type  of  a 
source  used.  The  greatest  stagnation  parameters  Po~14T07Pa  and  7o=9000K  with  an 
operation  time  of  10_3s  are  attained  in  the  case  of  employing  an  impulse  wind  tunnel 
plenum  chamber  as  a  hot  gas  source  [8],  while  7o~6800K  and  Po~108Pa  are  obtained 
when  application  is  given  to  a  region  behind  a  reflected  shock  wave  in  a  shock  wind 
tunnel  [9,10].  Various  arc  heater  designs  are  used  to  achieve  stationary  acceleration 
regimes,  for  example,  heaters  with  a  coaxial  electrode  arrangement  [11-13],  with  an 
extended  arc  [14]  etc.  The  temperature  and  pressure  ranges  in  them  are  Po=(2-20)-105Pa 
and  7o~3500-6000K,  the  heater  power  is  0.2-2MW  and  mass  flow  rate  through  the 
heater  is  0.02-0.25kg/s.  Nitrogen  and  air  are  used  as  a  test  gas.  The  stagnation 
parameters  in  stationary  sources  of  conducting  gases  can  be  increased  up  to  Po~2-107Pa, 
7o~5000K.  Some  experience  in  operation  with  similar  heaters  is  gained  both  in  Russia 
and  in  the  USA. 

To  provide  an  optimal  MHD-accelerator  operation,  the  accelerator  efficiency  defined  as 
the  ratio  of  additional  mechanical  energy  to  the  Joule  dissipation  function  r\=jBU/jE 
should  be  maximum.  The  efficiency  r\  increases  with  <jB2L,  therefore,  minimum 
possible  altitudes  should  be  simulated  using  a  maximum  possible  value  of  gB2L  when 
conduction  is  sufficient  to  have  a  satisfactory  current  density  in  the  acceptable  range  of 
electric  field  intensities.  Almost  in  all  publications  the  value  a=100mo/m  is  set  as  a 
lower  limit.  In  order  to  ensure  this  limit  it  is  necessary  to  inject  alkali  seeds  into  the  gas. 

2.2.  Seed  Injection  Methods 

The  method  of  seed  injection  and  provision  for  seed  uniformity  in  a  flow  govern,  to  a 
large  extent,  the  characteristics  of  a  MHD-accelerator.  There  are  several  ways  of 
injecting  seeds:  using  a  powder-type  K2CO1  gas  or  a  liquid  solution  of  K2CO3  or  KOH, 
applying  NaK  eutectics  and  supplying  K,  Cs.  Seeds  in  the  form  of  aerosols  are 
employed  for  impulse  systems.  The  fraction  of  total  seed  mass  is  varied  from  3%  to 
0. 5-0.3%.  The  tests  in  [1 1]  and  [12]  show  that  seeds  in  the  form  of  1  %NaK,  1  %K,  KOH 
do  not  exert  significant  effect  on  the  flow  characteristics.  In  a  number  of  previous 
publications,  erosion  properties  of  seeds  are  pointed  out  which  seem  to  be  caused  by 
incomplete  seed  evaporation.  Supplying  seeds  in  the  form  of  liquid  metal  aerosol 
through  atomizers  into  a  high-temperature  oxidizing  medium,  which  a  hot  air  is,  results, 
during  contacts  with  it,  in  formation  of  oxides  choking  the  duct.  To  avoid  it,  the  ducts 
should  be  blown  through  by  a  neutral  gas  prior  to  and  post  metal  supplying.  The  seed 
supply  area  should  be  chosen  so  that  a  uniform  seed  distribution  in  the  flow  be  provided 
and  seed  ingress  into  the  heater  discharge  zone  be  excluded.  To  do  this,  it  is  preferential 
to  use  a  special  mixing  chamber  [15]  with  dimensions  sufficient  for  metal  aerosol 
evaporation.  The  characteristic  metal  aerosol  evaporation  time  is  ~l(T3-10"2s.  In  future 
MHD-accelerators  it  seems  to  be  reasonable  to  inject  seeds  in  the  form  of  alkali  vapor. 
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Note  that  methods  of  providing  non-equilibrium  conduction  using  various  preionizers 
are  also  considered  in  the  literature.  Possible  methods  can  involve  photoionization, 
high-frequency  discharge,  radiometric  method,  glow  discharge  and  ionization  by 
electron  beam.  But  in  spite  of  a  30-year  discussion,  these  methods  did  not  yet  achieve 
the  level  of  practical  implementation  [16-18]. 

There  are  two  aspects  of  the  seed  effect  on  the  characteristics  of  aerodynamic  flows. 
Seed  distribution  in  a  gas,  seed  effect  on  thermodynamic  medium  properties  and 
especially  on  chemical  process  rates,  on  deactivation  of  vibrational  degrees  of  freedom 
of  nitrogen  and  oxygen.  A  uniform  seed  distribution  is  ensured  by  seed  supplying  into 
the  mixing  chamber  from  a  metering  device  through  atomizers  providing  its  fine 
spraying  (-10-2 0pm)  and  evaporation  before  entering  the  primary  nozzle.  The  NaK 
atom  distribution  uniformity  in  the  flow  can  be  evaluated  using  intensities  of  their 
resonance  fluorescence  excited  by  laser  radiation  with  a  resonance  wavelength. 
According  to  tests  carried  out  at  TsAGI  (Fig.4),  a  uniform  seed  distribution  in  a  flow  is 
quite  possible.  Because  of  a  low  ionization  potential  and  a  small  partial  seed  pressure  in 
thermodynamic  equilibrium  the  seed  presence  can  be  manifested  in  the  form  of  some 
additional  energy  sink  into  internal  degrees  of  freedom  within  the  range  of  temperatures 
and  pressures  at  which  physico-chemical  processes  in  air  (excited  vibrations  and  disso¬ 
ciation)  proceed  weakly.  The  seed  effect  on  entropy,  dissociation  degree,  sound 
velocity,  viscosity  is  insignificant,  while  specific  heat  capacity  and  heat  transfer 
coefficient  are  affected  to  a  greater  extent  [19].  The  seed  influence  on  deactivation  of 
vibrational  degrees  of  freedom  and  chemical  reactions  will  be  analyzed  below. 

2.3.  Primary  Nozzle 

The  investigation  of  primary  nozzle  characteristics  shows  that  the  most  optimal  range  of 
Mach  numbers  at  the  MHD-channel  entry  is  within  M=l. 8-2.2.  At  lower  Mach 
numbers,  there  occurs  some  flow  instability  at  the  MHD-channel  entry,  while  at  higher 
Mach  numbers  the  static  temperature  becomes  insufficient  to  provide  conduction,  and  a 
required  increase  in  the  stagnation  temperature  in  the  conducting-gas  source  results  in 
flow  contamination  by  the  electrode  materials  because  of  electrode  erosion. 

According  to  the  theoretical  and  experimental  data,  the  primary  nozzle  flow  in  all 
accelerators  subjected  to  study  is  in  thermodynamic  equilibrium.  In  [15],  use  is  made  of 
contoured  primary  nozzles  with  a  gradientless  uniform  velocity  field  at  the  exit. 

2.4.  MHD-Channel 

The  gas  flow  in  the  Faraday-type  accelerator  MHD-channel  differs  fundamentally  from 
flows  in  the  MHD-generator.  The  matter  is  that  the  MHD-channel  is  much  more 
stressed  in  terms  of  energy.  The  flow  at  the  channel  entry  features  a  supersonic  Mach 
number  exceeding  typical  values  for  the  generator  and  much  higher  stagnation  enthalpy 
(~2-3  times).  The  current  density  in  the  flow  is  respectively  higher  (/~50A/cm2)  at 
magnetic  induction  values  of  5=1. 5-5. 0T  and  more,  the  Hall  parameter  being  wexe~ 3- 
20.  Combining  the  flow  acceleration  and  temperature  rise  in  the  accelerator  results  in  a 
great  increase  in  the  stagnation  enthalpy  as  compared  to  the  initial  one  (4-6  times).  At 
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the  MHD-channel  exit,  the  Mach  flow  number  can  exceed  M>4.  Under  these 
conditions,  a  maximum  flow  temperature  is  attained  in  the  boundary  layer  near  the  wall 
(recovery  temperature).  Because  of  the  current  concentration  near  the  electrodes,  the 
Joule  gas  heating  at  low  temperatures  immediately  near  the  wall  brings  additionally  a 
high-temperature  and,  accordingly,  a  high-conducting  region  closer  to  the  surface.  The 
gas  density  in  this  region  is  much  less  than  that  in  any  other  flow  parts.  Because  the 
accelerating  Lorentz  force  is  a  volume  force  the  gas  in  this  region  acquires  a  greater 
velocity  than  the  main  flow  which  is  favorable  for  the  flow  pattern.  On  the  other  hand, 
the  same  phenomenon  increases  considerably  the  boundary  layer  displacement 
thickness  as  compared  to  an  ordinary  gas  dynamic  increase.  A  great  potential  drop 
across  the  cold  sublayer  of  the  boundary  layer  gives  rise  to  phenomena  similar  to 
electrical  break-down  followed  by  the  formation  of  microarcs.  Although  this 
phenomenon  is  well  known  in  the  extensive  literature  devoted  to  microarc 
investigations  in  MHD-generators  [20-21]  and  in  MHD-accelerators,  it  has  its  peculiar 
features  [22-24]. 

In  the  MHD-accelerator,  the  arcs  on  the  cathode  occur  and  concentrate  basically  at  its 
upstream  end.  Under  the  action  of  gas  dynamic  forces  they  are  carried  away  at  a  mean 
velocity  of  ~10m/s,  extinguish  and  then  occur  again  at  the  leading  edge.  A  characteristic 
current  in  an  individual  arc  is  5-1 0A.  The  quantity  of  microarcs  increases  with  rising 
the  current  strength  per  electrode.  The  anode  microarcs  are  stationary  to  be  placed  at  the 
joint  of  the  anode  trailing  edge  with  the  interelectrode  insulator  thus  causing  intense 
anode  erosion.  The  current  layers  of  high  conduction  lie,  as  a  rule,  behind  the  anode  arcs 
to  be  a  peculiar  kind  of  gas  electrodes.  This  phenomenon  is  similar  to  that  taking  place 
in  electric  discharge  glow  in  a  supersonic  flow  [25].  The  current  can  move  along  these 
layers  by  large  distances  thus  shunting  individual  electrodes  [24].  From  the  above 
outlined  phenomena  it  is  inferred  that  a  maximum  fine  MHD-channel  sectionalization  is 
desirable  and  the  microarc  discharge  motion  control  in  the  channel  is  obligatory.  The 
volt-ampere  characteristics  of  the  channel  discharge  gap  are  given  in  [26]  and  [27], 
They  are  not  volt-ampere  characteristics  in  the  classical  understanding  because  gas 
dynamic  flow  parameters  U,  T,  a,  P  vary  as  the  current  strength  rises.  Note  that  the  volt- 
ampere  characteristics  of  the  accelerator  obtained  at  AEDC,  USA  differ  qualitatively 
from  those  found  at  TsAGI,  Russia.  At  AEDC,  the  characteristics  feature  an  increasing 
behavior,  while  at  TsAGI  they  are  almost  flattened,  i.e.,  voltage  increases  slightly  with 
rising  the  current  strength.  The  lack  of  gas  dynamic  data  for  the  AEDC  accelerator 
makes  it  impossible  to  clarify  causes  of  the  difference. 

2.5.  MHD-Channel  Designs 

MHD-channel  designs  developed  at  different  organizations  differ  noticeably  to  depend 
on  a  chosen  type  of  the  conducting-gas  source  and  the  acceleration  process  duration. 
The  common  point  for  all  channels  is  the  Faraday  scheme  consisting  of  a  quadrangular 
channel  with  two  insulated  walls  and  two  electrode  walls  having  sectionalized 
electrodes  and  a  stand-alone  electric  power  source.  The  quantity  of  electrodes  in  the 
channel  is  from  3-5  to  120  at  different  organizations.  Extensive  studies  were  carried  out 
to  choose  materials  for  electrodes  and  insulators.  Most  suitable  materials  for  electrodes 
under  the  above  conditions  proved  to  be  copper  and  its  alloys,  while  boron,  silicon  and 
aluminium  nitrides  and  beryllium  oxide  deposited  on  the  cooled  surfaces  were  chosen 
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for  insulators.  In  stationary  facilities,  the  tests  were  conducted  at  static  pressures  of 
0.05-0.5 -105Pa  and  magnetic  induction  of  5=1. 5-2. 5T,  i.e.,  at  values  provided  by  the 

,  9  9  _ 1  _ O 

magnets  with  iron  core.  In  impulse  facilities  (t~10  -10  s),  the  pressures  were 
P~5.5-105Pa  and  5=4-5 T.  In  all  channels,  the  conductivities  were  within  the  range  of 
a~100-400mo/m,  current  densities  through  electrodes  for  stationary  and  impulse  MHD- 
accelerators  were  20-60A/cm2.  At y >60 A/cm2,  intense  electrode  erosion  was  noted.  It  is 
usual  to  point  out  a  mean  current  density  based  on  the  electrode  area  without 
consideration  of  the  current  concentration  at  the  electrode-insulator  joint.  The  electrodes 
and  insulating  walls  were  made  in  two  modifications:  with  water  cooling  and  without  it. 
In  the  modification  without  water  cooling  they  provide  heat  accumulation  in  their 
structures,  their  operation  time  being  restricted  by  attaining  the  insulator  surface 
temperature  corresponding  to  melting  or  evaporation.  In  publication  [15],  this  time  for 
uncooled  structure  is  -Is.  If  the  structure  is  cooled  by  water,  as  is  the  case  in  [11-12], 
the  time  is  ~10s  and  more.  In  similar  structures,  the  performance  of  cooled  copper 
electrodes,  as  well  as  of  cooled  insulators  made  from  a  thin  film  of  beryllium  oxide 
deposited  on  the  metal  is  good.  It  is  known  that  among  insulators  the  beryllium  oxide 
features  the  highest  heat  conduction  and  the  smallest  increase  in  conduction  as  the 
temperature  rises.  The  channels  were  equipped  with  probes  to  measure  pressure 
distributions  along  the  anode  and  cathode  walls,  as  well  as  transducers  for  voltage  and 
current  strength.  There  is  no  data  so  far  about  potential  measurements  inside  the  MHD- 
channel  of  the  accelerator.  In  almost  all  structures,  the  channel  is  made  slightly 
divergent  to  compensate  for  the  boundary  layer  growth,  the  range  of  divergence  angles 
being  0.3-1. 5°.  There  exist  channels  with  both  divergent  electrode  wall  only  and  all 
divergent  walls. 

2.6.  Electric  Power  Supply  for  Electrodes  and  Magnets 

50-100  stand-alone  electric  power  sources  providing  the  current  of  -50-100A  and 
voltage  -300-1 000 V  are  required  for  MHD- accelerator  electrodes,  the  insulators  being 
designed  for  voltage  of  l-5kV.  In  stationary  systems,  use  was  made  of  accumulator 
batteries,  high-voltage  transformers  with  disconnected  coil  winding  at  the  high-voltage 
side  in  accordance  with  voltages  required  for  individual  coils  connected  with  rectifying 
bridges,  individual  transformers  with  rectifying  devices  for  each  electrode  and  feedback 
systems  to  maintain  the  current  during  tests.  Capacitor  banks  were  used  for  electrodes 
and  magnets  in  impulse  systems.  The  electromagnets  in  stationary  MHD-accelerators 
were  energized  by  electric  generators.  At  AEDC,  a  project  was  developed  to  energize 
the  MHD-accelerator  from  the  MHD-generator,  but  this  project  was  not  realized  yet. 

The  table  given  below  summarized  the  main  characteristics  of  well-known  stationary 
MHD-accelerators  of  the  USA,  France  and  Russia. 
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Table  1.  MHD-accelerators  with  stationary  flows  for  wind  tunnels 


Organization 

Langley,  USA 

AEDC,  USA 

ONERA, 

France 

TsAGI,  Russia 

TsAGI, 

Russia 

Project 

1 

2 

3 

4 

5 

6 

n 

24 

117 

31 

L,  m 

0,3048 

0.77 

0.50 

2.54x2.54 

2. 5x2.5 

2.5x3. 5 

2.5x1.15 

| 

0.0068 

0.1 

0.04-0.25 

0.01-0.02 

1.14 

1.55-2.05 

1.4 

"2.5 

12 

16-20 

25-60 

25-30 

30-60 

30-40 

0.05 

0.5 

0.5 

0.2 

4400-5400 

2800-3200 

3000 

3000 

E9 

3-6 

4.5 

3.0 

200 

To,  K 

4000-6000 

3500 

3700 

4700 

Mb 

1.6 

1.6 

2.0 

2.0 

Uout,  m/s 

6000 

3800 

3200 

6000 

Gas 

Nitrogen 

Nitogen 

Air 

Air 

Seed 

l%Cs 

l%NaK 

7%KOH 

l%NaK 

(C2H5OH) 

■■ 

N,  kW,  heater 

55 

600 

200 

Nacc,  kW 

200 

1000 

2000 

n  -  number  of  electrodes;  L  -  accelerator  length;  F-  section  area;  Q  -  gas  mass  flow;  B 
-  magnetic  field  strength;  Po,  To  -  stagnation  parameters  in  the  conducting-gas  source; 
M  -  Mach  number  at  the  accelerator  entry;  U  -  velocity;  N  -  power. 
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3.  NUMERICAL  SIMULATION 

3.1.  Calculation  of  Gas  Flows  In  MHD-Channel 

The  gas  acceleration  process  in  the  MHD-channel  can  be  considered  from  two 
standpoints:  microscopic  and  macroscopic  ones.  The  microscopic  approach  analyzes  the 
electron  and  ion  behaviour  in  mutually  perpendicular  E  and  H  fields.  In  both  cases,  the 
accelerating  force  acting  on  the  gas  is  produced  by  ions  which  are  accelerated  by 
external  fields  within  the  interval  between  collisions  to  transfer  their  momentum  to 
neutral  particles. 

The  flow  in  the  Faraday-type  MHD-channel  with  sectionalized  electrodes  is  known  to 
be  an  essentially  three-dimensional  one.  Its  reduction  to  a  one-dimensional  canonical 
flow  derived  by  averaging  the  equations  of  magnetic  hydrodynamics  necessitates  strict 
assumptions.  Most  significant  among  them,  namely  the  absence  of  the  Hall  current 
along  the  channel,  current  non-uniformity  in  the  interelectrode  gap  and,  accordingly, 
increased  resistance  can  be  taken  into  account  by  introducing  effective  conductivity  <ieff 
depending  on  the  Hall  parameter  and  the  channel  sectionalization  degree  [28].  Voltage 
drop  values  in  the  boundary  and  electrode  layers  can  be  represented  in  the  form  of 
correction  coefficients  and  volt  equivalents  for  the  anode  and  the  cathode,  respectively 
[27].  In  this  case,  the  equations  for  gas  state  required  to  close  the  set  of 
magnetogasdynamic  equations  can  be  expressed  both  for  a  gas  in  thermodynamic 
equilibrium  in  the  form  of  H=H(P,T),  S=S(P,T)  and  for  a  thermodynamically  non¬ 
equilibrium  model  in  the  form  of  a  set  of  equations  of  physico-chemical  kinetics.  This  is 
also  the  case  for  transfer  coefficients  p,  k  etc. 

A  main  problem  of  experimental  investigations  lies  in  the  verification  of  all  above 
outlined  assumptions  and,  in  particular,  electrode  losses  and  boundary  layer  structure. 

A  set  of  magnetic  hydrodynamics  equations  with  the  above  assumptions  was  applied  at 
TsAGI,  Russia  to  analyze  the  gas  flow  in  the  MHD-accelerator  channel.  As  initial  and 
boundary  conditions  use  was  made  of  the  following  experimental  data:  gas  dynamic 
parameters  at  the  MHD-channel  entry  P,  T,  U,  M  and  distributions  along  the  channel 
P=P(x),  J=J(x),  U=U(x).  The  gas  dynamic  parameters  at  its  exit  obtained  experimentally 
served  for  verification  of  the  obtained  solution  (Fig.3).  This  procedure  is  described  in 
detail  in  [29-30],  Further  improvement  of  the  one-dimensional  calculation  technique  for 
MHD-channels  suggests  the  consideration  of  non-equilibrium  in  the  MHD-channel  and 
a  more  detailed  analysis  of  energy  transport  channels  from  the  electric  field  to  the 
kinetic  flow  energy  assuming  that  the  electrons  accelerated  in  the  electric  field  transfer 
their  energy  to  vibrational  degrees  of  freedom  of  air  molecules  and  seed  atoms. 
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Test  results: 


Analytic  results: 


static  pressure  P-1 0‘5Pa 
electrode  voltage  V-10‘3B 
current  through  electrodes  J-10'2A 

velocity  defined  by  measuring  the  induced  voltage  value  U-lO^m/s 


pressure  P-10~5Pa 
temperature  T-10'4K 
velocity  U-10‘4m/s 
Mach  number  M  - 1 0'1 


Light  points  show  the  MHD-accelerator  operating,  dark  points  show  Solid  line  =  the  accelerator  operating, 
the  absence  of  MHD-acceleration  dash  line  =  without  MHD-acceleration. 


Fig-3.  Distribution  of  electric  and  gas-dynamic  parameters  along  the  accelerator 

channel 


The  schematic  of  energy  transport  from  the  electric  field  to  internal  degrees  of  freedom 
is  given  below. 
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The  losses  of  applied  energy  for  radiation  and  for  heat  fluxes  towards  the  electrodes  and 
the  walls  are  accounted  for  through  the  loss  coefficient  obtained  experimentally.  The 
calculations  performed  show  that  the  vibrational  temperature  in  the  TsAGI  MHD- 
channel  proves  to  be  sufficiently  close  to  the  translational  one  because  of  increasing  the 
deactivation  rate  of  the  vibrational  temperature  due  to  O  and  N  atoms  in  the  flow. 
Comparison  of  calculated  output  parameters  for  the  same  experimental  distribution  of 
currents  J(x),  voltages  between  the  electrodes  U(x)  and  static  pressures  P(x)  along  the 
MHD-channel  reveals  that  the  differences  in  the  flow  velocity  and  temperature 
calculated  under  assumption  of  thermodynamic  equilibrium  and  without  it  are  not  very 
great,  while  the  difference  in  the  chemical  gas  composition  is  appreciable.  Particular 
differences  for  the  TsAGI  MHD-accelerator  are  presented  in  [15].  It  is  significant  that 
as  the  static  pressure  in  the  MHD-channel  rises  the  degree  of  deviation  from  the 
equilibrium  decreases.  However,  as  the  pressure  rises,  the  heat  flux  q  towards  the 
electrode  and  the  insulating  walls  of  the  MHD-channel  increases  proportionally  [31]. 

The  papers  [32-33]  contain  the  experimental  data  on  convective  heat  flux  components 
and  the  heat  flux  part  caused  by  the  current.  The  convective  component  is  determined 
using  heat  flux  values  in  upstream  and  downstream  discharge  zones,  respectively,  for 
the  anode  and  cathode  walls.  The  flow  parameters  (p,  u)  required  for  the  calculations 
are  obtained  under  the  assumption  of  thermodynamic  equilibrium.  It  is  found  that  the 
convective  heat  flux  component  is  satisfactorily  described  by  the  Spolding-Chee  theory 
for  a  turbulent  boundary  layer  [32],  though,  according  to  Reynolds  numbers  attained  in 
the  MHD-channel,  the  flow  lies  in  the  transitional  region  between  laminar  and  turbulent 
flows.  The  boundary  layer  disturbances  due  to  microarc  discharges  seem  to  cause  fast 
boundary  layer  turbulization. 

An  intense  local  gas  heating  in  the  microarc  at  y>500A/cm2  can  cause  a  boundary  layer 
separation  followed  by  formation  of  weak  shock  waves.  Some  regions  with  a  spatial 
discharge  division  and  with  a  high  conduction  occur  in  the  shock  waves  along  which  the 
current  flows.  The  appearance  of  local  boundary  layer  separation  regions  leads  to  an 
increase  in  the  aerodynamic  drag  of  the  MHD-channel.  Based  on  the  experimental  data 
obtained  at  TsAGI,  the  friction  coefficient  can  be  taken  close  to  C/~5-lCf3. 

It  is  also  necessary  to  account  for  possible  three-dimensional  effects  caused  by  the 
development  of  secondary  flows  in  extended  MHD-channels,  the  influence  of  electrode 
shunting  by  anode  "gas"  electrodes,  various  effects  of  thermodynamic  non-equilibrium 
before  we  can  state  that  the  theory  of  MHD-accelerator  is  constructed. 

3.2.  Flow  Calculation  in  the  Secondary  Nozzle 

The  MHD-accelerator  exit  is  equipped  with  secondary  nozzles.  The  gas  expansion  level 
in  them  is  specified  by  test  objectives  and  capabilities  of  the  pumping  system.  There  are 
almost  no  information  about  gas  flows  in  secondary  nozzles  in  previously  mentioned 
publications  of  the  USA  and  France.  In  investigations  carried  out  at  TsAGI,  Russia  the 
gas  expansion  level  in  the  secondary  nozzle  was  FOM/F„CC>100.  The  investigations 
involved  pressure  distributions  Pst  along  the  nozzle  and  total  pressure  field  behind  the 
normal  shock  at  its  exit  Pq  .  Direct  application  of  experimental  data  to  determine  gas 
dynamic  parameters  presents  certain  difficulties  because  of  uncertainty  in  the 
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thermodynamic  flow  model.  It  was  suggested  in  previous  papers  [29]  that  the  flow 
corresponded  to  a  frozen  one  with  an  adiabatic  exponent  y  =  5+2^£,(-/3  +  2^£,/ 
determined  for  the  mixture  of  uniatomic  and  biatomic  gases,  where  is  the  mole 
fractions  of  biatomic  02,  N2  and  NO  molecules.  Their  values  were  chosen  for 
thermodynamically  equlibrium  gas  flows  from  the  MHD-channel.  A  more  detailed 
analysis  made  in  [34]  with  due  consideration  of  two-temperature  dependence  of  rates  of 
recombination  and  deactivation  of  vibrational  degrees  of  freedom  shows  that  the  seed 
accelerates  the  deactivation  of  vibrational  degrees  of  freedom  of  N2  resulting  in  some 
increase  in  physical  gas  velocity,  growth  in  the  static  temperature  and  decrease  in  the 
Mach  number.  The  experimental  Mach  number  at  the  secondary  nozzle  exit  is  M~15  at 
a  flow  velocity  of  F~7800m/s  and  temperature  of  -500K,  P^=20-50Pa.  Measurements 
of  distributions  of  Pst  along  the  nozzle  and  fields  of  Pq  make  it  possible  to  verify  the 
validity  of  the  thermodynamic  flow  model  used  in  the  calculations.  The  uniformity  level 
for  fields  of  Pq  at  the  nozzle  exit  depends  on  flow  characteristics  in  the  MHD-channel, 
in  particular,  on  the  Hall  current.  When  its  value  is  great  the  pressures  on  the  cathode 
and  anode  walls  differ  considerably,  while  the  flow  field  becomes  non-uniform. 
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4.  EXPERIMENTS  WITH  MHD  WT 

4.1.  Experimental  Investigation  of  Flow  Parameters 

The  aerodynamic  investigations  require  reliable  gas  dynamic  and  physico-chemical 
flow  parameters,  as  well  as  respective  means  of  flow  visualization  and  measurements  of 
forces,  pressures  on  models  and  heat  fluxes  towards  them  etc.  A  main  governing 
parameter  is  the  flow  velocity  at  the  accelerator  and  secondary  nozzle  exits.  There  exist 
several  methods  of  its  measurement.  A  simple  and  low-cost  way  implies  measuring  the 
potential  induced  on  one  of  the  last  non-energized  electrode  pair  U=uBh,  where  h  is  the 
distance  between  the  electrodes,  u  is  the  flow  velocity.  It  is  a  common  practice  at 
TsAGI,  Russia  to  use  the  gas  dynamic  method  based  on  comparison  of  two 
measurements  of  Pq  at  one  and  the  same  point,  i.e.,  first  measurement  when  the  velocity 
value  can  be  derived  from  simple  gas  dynamic  relationships  under  assumption  of 
adiabatic  and  isentropic  flow  nature  at  known  values  of  Po  and  To,  and  the  second 
measurement  when  the  gas  flow  is  subjected  to  some  action  and  Joule  heat  is  supplied 
to  it.  In  this  case,  assuming  that  the  profiles  of  pressures  Pq  are  similar  and  gas  mass 
flow  rates  in  the  channel  with  and  without  MHD-acceleration  are  constant  we  obtain: 


pi 

£02  with  MHD 
p 

1  01  without  MD 


FP2U2  1' 


8  / with 


MHD 


A 


PpiMj 


8  f without  MBD 


,«2 

U\ 


Another  method  developed  for  this  purpose  at  TsAGI  is  the  method  of  measuring  the 
velocity  based  on  the  Doppler  effect.  In  this  method,  K  or  Na  seed  atom  glow  is  excited 
by  laser  radiation  with  a  resonance  frequency  at  a  specified  flow  point,  a  wavelength 
shift  being  measured  by  a  Fabry-Perot  interferometer.  The  accuracy  of  measuring  the 
velocity  of  ~6000m/s  is  ~3-5%.  Some  versions  of  this  device  were  also  developed  to 
obtain  panoramic  flow  velocity  values  at  its  specified  section,  i.e.,  to  obtain  a  velocity 
field  in  any  test  section  part  or  at  the  nozzle  exit. 

The  paper  [35]  considers  the  results  of  measurements  of  flow  velocity  fields  at  the 
channel  exit  by  analyzing  the  motion  of  an  artificial  or  a  natural  non-uniformity,  but  the 
accuracy  of  this  method  is  small.  The  velocity  value  can  be  calculated  using  a  certain 
thermodynamic  flow  model  and  applying  experimental  distributions  of  pressures, 
current  strengths  and  voltages  along  the  channel  and  specified  parameters  at  its  entry. 
Agreement  of  velocity  values  obtained  by  different  methods  used  at  TsAGI  proves  to  be 
good. 

It  is  impossible  to  determine  the  chemical  gas  composition  at  the  nozzle  exit  by  direct 
methods.  Among  indirect  methods,  the  simplest  method  is  determining  y=y(^)  which  is 
possible  using  static  pressure  distributions  along  a  hypersonic  nozzle  for  known 
pressure  field  of  Pq  and,  accordingly,  displacement  thickness.  The  vibrational 
temperature  level  for  N2,  02  and  NO  molecules  can  be  found  by  methods  of  emission 
spectroscopy  with  an  acceptable  accuracy  [36].  A  more  detailed  information  about 
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thermodynamic  flow  state  including  level-by-level  kinetics  can  be  obtained  using 
CARC  methods,  but  as  far  as  is  known,  this  method  was  not  applied  to  analyze  flows  in 
MHD-accelerators . 

4.2.  Application  of  Hypersonic  MHD-Gas  Acceleration  Wind  Tunnels  to 
Aerodynamic  Experiments 

Unfortunately,  the  developments  of  hypersonic  MHD-gas  acceleration  wind  tunnels 
were  ceased  in  the  USA  and  France  at  the  initial  stage  and  no  results  of  their  application 
to  solve  problems  of  aerodynamics  were  published. 

Visualization 

At  TsAGI,  similar  investigations  were  started  in  the  80s.  They  included,  first  of  all, 
developments  of  test  procedures  as  applied  to  visualization  of  flow  patterns  over  space 
vehicle  models  during  their  atmospheric  reentry  and  investigations  of  pressure  and  heat 
flux  distributions  on  vehicle  surfaces,  theoretical  investigations  of  the  flow  pattern  in 
the  presence  of  physico-chemical  reactions  in  the  boundary  layer  near  models  and  at 
considerable  "frozen"  dissociation  level  and  vibrational  temperature  of  molecules  in  a 
ffee-stream  flow. 

The  application  of  Na,  K  seeds  in  the  flow  offers  radically  new  possibilities  based  on 
using  the  effects  of  resonant  fluorescence  and  abnormal  dispersion  to  visualize  flows 
over  hypersonic  vehicle  models.  Employing  laser  radiation  with  a  resonance  wavelength 
of  X=5890A  in  the  laser  knife  method  and  optical  light  filters  which  exclude  intense 
natural  gas  radiation  enables  the  shock  wave  position  and  shape,  internal  flow  pattern 
vortices,  flow  separation  and  attachment  zones  to  be  obtained  (Fig.4).  As  a  result  of  the 
application  of  the  resonance  laser  radiation  as  a  light  source  in  the  shadow  device  and 
the  abnormal  dispersion  effect,  the  device  sensitivity  increases  by  a  factor  of  100  and 
more  and  intense  natural  gas  radiation  is  eliminated  [37]. 

Seed  Effects 

A  detailed  analysis  of  the  seed  influence  on  the  flow  over  a  model  shows  that  the  seed 
does  not  exert  any  effect  on  the  thermodynamic  gas  state  in  the  critical  flow  line,  but  it 
changes  the  position  of  weak  shock  waves  at  great  distances  from  the  model  due  to 
accelerated  deactivation  of  vibrational  degrees  of  freedom  of  N2  and  O2. 

The  primary  seed  effect  consists  in  increased  electron  concentration  and  radiation 
intensity  in  resonance  lines.  Heat  flux  values  are,  however,  much  smaller  than 
convective  ones,  therefore,  they  do  not  essentially  influence  the  tests  data. 

Code  Validation 

Owing  to  well-studied  techniques  of  determining  gas  dynamic  flow  parameters  and  flow 
visualization  methods  it  became  possible  to  verify  calculation  codes  and 
thermodynamic  models  in  hypersonic  gas  flows,  as  well  as  to  conduct  comparative  tests 
in  different  wind  tunnels  to  evaluate  real  gas  effects.  The  experiments  yielded  pressure 
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distributions  on  surfaces  of  simple  models  and  the  data  on  shock  wave  position  and 
shape  over  them.  It  was  revealed  that  a  gas  glow  region  fixed  visually  near  the  model, 
which  was  taken,  as  a  rule,  as  a  shock  wave,  was  considerably  larger,  almost  2  times 
than  its  true  position  obtained  using  the  shadow  method  including  the  abnormal 
dispersion  effect.  This  effect  can  be  explained  by  the  displacement  of  radiation  and  the 
diffusion  of  high-energy  electrons  from  a  high-temperature  region  (7M0000K)  behind 
the  shock  wave  further  downstream.  This  phenomenon  was  known  earlier  only  in 
investigating  shock  waves  caused  by  nuclear  explosion  and  meteorites. 


Fig.4.  Normal  cross  section  of  shock  wave  at  Bor-4  aircraft  visualizated  by  resonance 

fluorescence  method. 


The  flow  visualization  region  behind  the  MHD-accelerator  is  also  of  great  interest  for 
fundamental  investigations  of  gas  properties  at  very  high  temperatures  and  pressures 
(r0>10000K,  P0>105Pa). 

Comparison  of  theoretical  and  experimental  data  on  pressure  distributions  and  shock 
wave  positions  and  shapes  made  it  possible  to  assess  the  validity  of  some  or  another 
flow  model,  as  well  as  the  influence  of  thermodynamic  medium  properties. 

Non-Perfect  Gas  Effects 

How  dissociation  and  vibrational  excitation  of  molecules  in  a  free-stream  flow,  as  well 
as  alkali  seed  disturb  the  flow  pattern  over  models  and  what  is  the  level  of 
correspondence  of  wind  tunnel  results  to  full-scale  data  are  the  problems  of  crucial 
importance  which  affect  further  developments  of  wind  tunnels  with  MHD-accelerators. 
To  study  these  problems,  gas  dynamic  and  physico-chemical  parameters  of  a  flow  over 
a  sphere  were  calculated  for  full-scale  and  wind  tunnel  conditions  assuming  the  equality 
of  total  stagnation  enthalpy,  dynamic  pressures  and  binar  similarity  parameters.  The 
calculations  were  carried  out  for  a  standard  facility  operation  regime  [29]  which 
corresponded,  according  to  the  above-outlined  conditions,  to  a  flight  at  an  altitude  of 
62.5km  and  at  a  velocity  of  7.430m/s.  The  complete  description  of  the  calculation 
technique  and  the  results  obtained  are  given  in  [38].  It  follows  from  the  results  obtained 
that  the  stand-off  distance  of  the  shock  wave  from  the  sphere  in  the  wind  tunnel  is 
slightly  greater  than  that  in  full-scale  conditions  due  to  the  presence  of  oxygen  atoms 
amd  a  higher  level  of  excitation  of  vibrational  degree  of  freedom  of  molecules  in  a  free- 
stream  flow  over  the  model.  The  vibrational  temperature  profiles  in  full-scale  and  wind 
tunnel  conditions  are  similar  for  N2  and  O2  molecules,  while  for  NO  the  difference  at  the 
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wave  front  is  considerable  to  become  similar  near  the  body  surface.  The  excitation  of 
vibrational  degrees  in  the  wind  tunnel  and  in  full-scale  conditions  is  adequate. 

The  analysis  of  mass  concentration  profiles  for  O,  N,  NO  components  across  the  shock 
layer  reveals  that  under  full-scale  conditions  the  chemical  reactions  proceed  slower 
because  of  vibrationally  dissociative  interaction. 

Comparison  of  electron  densities  shows  that  in  the  MHD  facility  the  electron  energy  at 
the  forward  critical  point  exceeds  full-scale  values  about  three  times,  while  in  the  vicinity 
of  the  side  surface  the  agreement  is  good.  The  calculation  data  on  aerodynamic  heating  of 
the  body  surface  demonstrate  a  satisfactory  agreement.  A  greater  effect  of  catalytical 
surface  properties  for  the  MHD-conditions  is  associated  with  the  presence  of  dissociated 
atoms  and  vibrationally  excited  molecules  in  a  free-stream  flow. 

Good  agreement  of  calculated  and  measured  positions  of  the  shock  wave  front  near  the 
body  confirms  rather  reliably  the  validity  of  the  above-considered  results.  It  is  inferred 
from  the  analysis  carried  out  that  an  increased  degree  of  dissociation  and  vibrational 
excitation  of  molecules  in  a  free-stream  flow  changes  in  a  certain  way  the  flow 
conditions.  In  view  of  this,  it  is  desirable  to  choose  the  facility  operation  conditions  so 
that  their  values  be  minimum.  Nevertheless,  the  facility  with  already  attained  values  is 
suitable  for  testing  industrial  models  in  a  hypervelocity  air  flow. 

The  flow  over  the  Bor-4  vehicle  model  was  investigated  most  comprehensively. 
Comparison  of  calculated  and  experimental  data  on  shock  wave  position  and  shape  near 
this  model  made  it  possible  to  assess  the  significance  of  non-equilibrium  effects  in 
flows  [40].  The  distributions  of  heat  fluxes  on  the  model  surface  for  the  space  vehicle 
Bor-4  were  in  good  agreement  with  respective  flight  data. 

Scramjet  Simulation 

Of  great  importance  in  investigations  carried  out  using  this  facility  was  the  reproduction 
of  gas  dynamic  parameters  corresponding  to  the  conditions  in  the  scramjet  combustion 
chamber  in  the  facility  test  section  [39,40].  The  paper  [40]  presents  the  distributions  of 
heat  fluxes  close  to  full-scale  values  g=5-10kW/cm2  towards  such  scramjet  combustion 
chamber  components  as,  for  example,  fuel  pylons,  and  evaluates  the  non-equilibrium 
influence  behind  the  shock  wave  ahead  of  the  pylons.  The  results  obtained  are  compared 
with  the  Fay-Riddel  theory  including  possible  corrections  for  viscosity  effects  according 
to  the  Chang  theory.  It  is  shown  that  the  flow  behind  the  shock  wave  is  in  thermodynamic 
non-equilibrium.  In  the  tests  carried  out,  Mach  numbers  in  the  combustion  chamber,  flow 
velocity  V,  temperature  Ts,  (but  the  pressure  level  was  less  almost  10  times)  were 
reproduced.  The  total  enthalpy  values  proved  to  be  higher  than  those  calculated  using  the 
gas  dynamic  parameters  because  of  a  considerable  oxygen  dissociation  level. 

In  the  70s  the  problems  of  atmospheric  vehicle  reentry  were  solved  at  that  stage  by  using 
much  simpler  and  low-cost  facilities:  hot-shot  wind  tunnels,  shock  wind  tunnels  with 
reflected  shock  waves,  with  heavy  or  light  pistons,  classical  wind  tunnels  with  the 
application  of  various  corrections  for  gas  reality,  viscosity  and  compressibility  effects  and 
by  employing  computational  fluid  dynamics  methods  to  much  more  extent.  Early  in  the 
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90s,  the  interest  to  MHD-accelerators  was  rekindled  because  of  the  development  of 
transatmospheric  vehicles  (TAV  vehicles).  If  tests  on  models  of  a  smaller  scale  (1:10-25) 
during  a  short  time  period  of  x~10s  are  sufficient  for  investigating  airframes,  input 
engine  diffusors  and  integrated  engine  nozzle-airframe  configuration,  i.e.,  the  available 
test  facilities  can  be  used,  the  optimization  of  the  TAV  engine  necessitates  full-scale  tests 
during  a  rather  long  time.  There  are  no  adequate  aerodynamic  test  rigs  for  this  purpose 
and  hypersonic  MHD-gas  acceleration  wind  tunnels  prove  to  be  just  the  only  facilities 
essentially  suitable  for  similar  investigations  [41-42].  But  up  to  now  a  number  of 
fundamental  problems  concerning  the  influence  of  seeds  and  deviations  of  the  chemical 
compositions  of  free-stream  gases  from  natural  air  on  physico-chemical  processes  and  on 
oxygen  burning  remain  unsolved.  It  follows  from  the  calculations  that  the  mole  fraction  of 
"O"  atoms  (radicals)  in  the  free-stream  flow  more  than  10%  changes  considerably  the 
hydrogen  burning  process,  induction  time,  equilibrium  fraction  of  water  vapors,  as  well  as 
parameters  P,  p,  T  along  the  combustion  chamber.  Therefore,  this  facility  is  of  little  use 
for  testing  scramjets  because  an  acceptable  deviation  must  be  less  than  5%.  The  presence 
of  NO  in  the  flow  to  the  extent  of  ~5-5%  and  of  Na,  K  atoms  in  amounts  of  1%  does  not 
essentially  influence  the  burning  process. 


Wind  tunnel 
V  “  6.14  km/s 
p  =  3.23-lO^kg/m3 
Rw  ~  0.02  m 


=  0.94  -104K 
T™  =  2.21  •  103K 
AT  =  0.38  •  103K 


Fig.5.  Tm  isolines  in  the  shock  layer 


Operational  Map 

The  paper  [42]  considers  the  parametric  calculations  of  MHD-accelerator  under  the 
assumption  of  a  thermodynamical  equilibrium  in  the  channel  to  find  optimal  parameters 
sufficient  to  solve  the  scramjet  parameters.  The  calculations  were  performed  at  constant 
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stagnation  parameters  Po~2T07Pa,  7b=4700K  and  gas  mass  flow  rate  of  ~20kg/s,  and 
variable  magnetic  induction  S,  current  density  j,  Mach  number  at  the  channel  entry  and 
channel  wall  angles.  The  secondary  nozzle  parameters  were  chosen  so  that  P«=5-104Pa, 
Tsf= 2642K,  t/=6008m/c  were  attained  at  its  exit,  i.e.,  the  values  at  the  chamber  entry  at 
M=20.  According  to  this  paper,  it  is  quite  possible  to  obtain  flow  parameters  close  to 
required  ones  using  the  MHD-accelerator  at  Z?=12-16T.  In  the  paper  [43],  the  MHD- 
accelerator  was  calculated  at  the  same  initial  parameters  based  on  the  TsAGI  data  for  a 
thermodynamically  non-equilibrium  gas  and  with  the  consideration  of  a  possible  increase 
in  the  electron  temperature  as  compared  to  the  gas  temperature.  According  to  the  data 
obtained,  the  non-equilibrium  phenomenon  is  favorable  in  this  case  for  flow  parameters 
and  the  dissociation  level  at  the  channel  exit  decreases  almost  2  times.  Comparison  of 
parameters  realized  in  the  scramjet  combustion  chamber  and  at  the  secondary  nozzle  exit 
reveals  that  the  above  stated  problem  can  be  solved  at  Mach  flight  numbers  M=  15-20. 
When  this  flow  expands  in  the  secondary  nozzle  up  to  static  temperatures  corresponding 
to  atmospheric  conditions  it  is  possible  to  attain  in  the  test  section  full-scale  Mach 
numbers  (M=20-25),  full-scale  flight  velocities  at  pressures  corresponding  to  flight 
altitudes  of  H-6 5-80km.  The  resulting  flow  differs  slightly  from  the  real  one  by  the 
presence  of  NO  molecules  (5-6%)  and  O  atoms  (2-2.5%),  but  such  percentage  does  not 
influence  the  flow  over  the  vehicle  (Fig.6). 
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Fig.6.  Regions  of  flight  trajectory  modelling  in  wind  tunnels  of  different  classes.  1 . 
Missile  start  trajectory,  2.  Atmospheric  entry  trajectory  for  ballistic  heads,  3.  Flight 
corridor  for  orbiter  with  scramjet,  a)  Continuous  wind  tunnel,  b)  Intermittent  wind 
tunnel,  c)  Arc  heating  wind  tunnel,  d)  Hot-shot  wind  tunnel 


Future  Facility 

Comparison  of  this  MHD-channel  with  that  used  at  TsAGI  shows  that  the  power  applied  to 
the  new  channel  should  be  increased  1000  times,  the  pressure  100  times,  the  magnetic 
induction  5  times  and  the  channel  length  4-5  times,  the  relative  MHD-channel  length 
(number  of  gauges)  decreases  by  a  factor  of  3-4.  The  construction  of  a  similar  MHD- 
channel  poses  a  number  of  essentially  new  problems,  the  main  of  them  is  the  development 
of  a  serviceable  structure  designed  for  convective  heat  fluxes  towards  its  walls  of  q= 20- 
60MW/m2,  to  which  should  be  added  heat  fluxes  due  to  the  electrode  phenomena  qe~  10- 
20MW/m2,  thus  providing  a  specified  velocity  field  uniformity  level  at  the  nozzle  exit  etc. 
The  MHD-channel  structure  should  be  cooled  to  operate  within  r>10-60s.  It  is  unlikely  to 
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develop  a  similar  structure  because  the  insulator  surface  temperature  during  the  start  attains 
the  melting  temperature  for  t~1(T3s,  while  the  copper  electrode  surface  temperature  for 
t~5-10_1s.  Thin  insulating  film  coatings  deposited  on  the  copper  are  also  not  effective.  The 
only  way  is  to  create  a  gas  screen  from  a  light  gas  (helium)  enabling  the  heat  flux  to  the 
walls  to  be  decreased  considerably.  There  are,  however,  no  data  on  the  MHD-operation 
with  a  gas  screen.  The  serviceability  of  the  electrode-insulator  joint  can  be  ensured  if  some 
means  of  controlling  the  microarc  motion  on  the  electrode  surface  are  developed  to  prevent 
microarc  concentrations  at  the  electrode  end.  The  uniformity  of  velocity  fields  and  of 
chemical  air  composition  at  the  secondary  nozzle  exit  dictates  the  elimination  of  the 
influence  of  the  Hall  current,  current  stubs  due  to  microarcs  and  respective  three- 
dimensional  effects.  The  development  of  remaining  structural  components  of  a  hypersonic 
wind  tunnel  (arc  heater,  seed  supply  system,  primary  and  secondary  nozzles,  test  section, 
super-  and  subsonic  diffusors,  cooler  and  exhaust  system)  does  not  presents  considerable 
difficulties  except  for  the  magnetic  system  designed  for  B=\ 2T.  Thus,  the  development  of  a 
hypersonic  MHD-acceleration  wind  tunnel  capable  of  reproducing  correctly  the  scramjet 
combustion  chamber  entry  conditions  and  flight  conditions  at  65-75km  is  an  extremely 
complicated  problem  to  be  comparable  with  the  problem  of  developing  the  scramjet  itself. 
For  this  purpose,  it  is  necessary  to  attract  specialists  of  a  very  wide  spectrum  and  to  apply 
the  achevements  which  have  gained  in  high-temperature  magnetic  gas  dynamics  and 
hydrodynamics  for  the  lattest  30  years. 
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CONCLUSION 


1.  Hypersonic  MHD-gas  acceleration  wind  tunnels  developed  by  the  present  time 
are  used  successfully  to  solve  practical  problems  of  hypersonic  aerodynamics,  to 
verify  computational  codes  and  calculated  thermodynamic  medium  models  for 
non-equilibrium  gas  flows. 

2.  In  principle,  it  is  quite  possible  to  construct  a  MHD-accelerator  to  test  scramjets 
for  TAV  vehicles.  This  problem  is  comparable  with  the  development  of  the 
scramjet  itself.  It  will  demand  much  time  and  considerable  funding. 
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